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ABSTRACT 


This  study  examines  the  formation  of  a  front  over  a 
mountain  due  to  a  pre-existing  deformation  field.   The 
dependent  variables  are  assumed  to  be  independent  of  x.   The 
hydrostatic  Boussinesq  equations  are  used  with  no  diffusion 
of  heat  or  momentum.   A  modified  coordinate  system  similar 
to  Phillips'  sigma  system  is  used.   The  model  is  bounded  at 
the  top  by  a  rigid  plane  and  periodic  boundary  conditions  are 
used  in  the  horizontal.   The  solutions  show  the  formation  of 
a  front  within  a  finite  period  of  time  that  tilts  toward  the 
cold  air.   The  frontogenesis  experiments  with  a  mountain  of 
wavelength  3600  km  showed  negligible  topographic  effects. 
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I.   INTRODUCTION 

Williams  and  Plotkin  (19  68)  have  shown  that  atmospheric 
frontogenesis  occurs  for  large  time  limits  using  the  quasi- 
geostrophic  equation.   Their  results  are  unrealistic  in  the 
respect  that  the  frontal  zone  does  not  tilt  with  height.   A 
comparison  of  the  linear  and  nonlinear  hydrostatic  primitive 
equations  by  Williams  (1972)  shows  that  the  linear  solutions 
are  close  to  the  solutions  of  the  quasi-geostrophic  equations. 
The  nonlinear  solutions  are  more  realistic  and  show  the  forma- 
tion of  frontal  discontinuities  in  a  finite  period  of  time. 
Hoskins  and  Bretherton  (19  72)  showed  that  dynamic  processes 
in  addition  to  horizontal  deformation  can  produce  a  front. 
Bannon  (1983)  and  (1984)  investigated  quasi-geostrophic  and 
semi-geostrophic  frontogenesis  over  topography.   Monotonically 
sloping  topography  beneath  a  deformation  leads  to  the  forma- 
tion of  a  surface  discontinuity  in  the  absence  of  an  initial 
horizontal  thermal  gradient.   A  front  associated  with  a 
translating  deformation  field  over  an  isolated  orographic 
feature  is  affected  by  the  topographically  forced  divergence 
field  as  the  front  passes  over  the  moutnain.   Convergence 
will  increase  frontogenesis  and  divergence  will  decrease  it. 

The  purpose  of  this  study  is  to  investigate  the  effect 
of  topography  on  a  front  which  arises  from  deformation. 
Periodic  boundary  conditions  are  introduced  because  they 
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simplify  the  mountain  flow  simulation.   The  deformation  is 
independent  of  height.   The  Boussinesq  equations  are  used 
in  which  the  dependent  variables  are  independent  of  x.   A 
coordinate  system  similar  to  that  of  Phillips  (1957)  is  used. 
The  equations  are  solved  numerically  and  output  is  shown  in  the 
form  of  vertical  cross-sections  in  the  y-C  plane. 
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II.   BASIC  EQUATIONS  '  '      " 

The  basic  equations  used  in  this  study  are  the  Boussinesq 
equations  which  neglect  the  compressibility  of  the  atmosphere. 
The  hydrostatic  approximation  is  made  and  heating  and  friction 
are  neglected.   These  equations  will  be  applied  to  a  model 
in  which  the  flow  is  bounded  above  by  a  rigid  horizontal 
plane. 

The  Boussinesq  equations  [Ogura  and  Phillips,  1962]  may 
be  expressed  in  the  z-coordinate  system  as 


ll  +  v-Vv  +  w  P-     =   -V$  -  f  kxv  ,  (1) 

d  t  d  Z 


11  ,-.ve  .wff  =  0  ,  (2) 


If  =  ej^  '  <^' 


Vv  +  1^  =   0  ,  (4) 

a  Z 


where 


K 

0   =   T(-^)   -  0_  ,  (5) 


0   p  Pq       0 


and 
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7T     =   -  ^  Z  +  C    .  (7) 

0      e^     p 


In  this  model,  the  Coriolis  parameter  is  constant,  and  the 
only  quantity  which  varies  with  longitude  is  $. 

To  simplify  the  finite  differencing  of  the  equations,  a 
modified  coordinate  system  is  introduced  in  which  the  lower 
boundary  surface  becomes  a  coordinate  surface.   This  system 
permits  the  application  of  the  lower  boundary  condition  at 
the  surface  of  the  mountain.   The  coordinate  system  is  similar 
to  the  sigma  system  introduced  by  Phillips  (1957). 

Following  Walton  (1968)  and  DeBoer  (1970)  ,  the  vertical 
coordinate,  zeta,  is  defined  as 


z  -  z 
?   =  —^r-^    ,  (8) 


where  D  and  z   are  defined  in  Fig.  1  and  Z  =  D  -z  .   With 

s  ^  s 

these  definitions,  i,   =   0   a.t   the  lower  boundary  and  c  =  1  at 
the  upper  boundary. 

Transformation  of  equations  (1)  through  (4)  from  the 
z-coordinate  system  into  the  ^-coordinate  system  yields 


|v  ^  -.^-  ^    |v   ^   _  ^^_  |i(izl)v,  z^  -  f  kx  V  ,  (9) 

dt.  dC  d(^Zns 


II  +  v-ve  +  ?  H  =  0  ,  (10) 


H  =  ^^«  '  <ii' 
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V-  (Zv)  +  Z  ll-   =   0  . 


(12) 


In  equations  (9)  through  (12) ,  V  is  the  horizontal  gradient 
operator  on  zeta  surfaces. 

It  is  desirable  to  put  the  equations  in  flux  form  in  order 
to  satisfy  certain  computational  properties.   Equations  (9) 
through  (12)  become 


3  (M  )  Pi   •  ;^6 

^.^   +  i-(M  U)  +  47-(?uZ)   =   -  Z4^  +  fM   , 
9t      8y   y      8^  ^  8x      y 


(13) 


d  (M  )  . 

^t^  ^  i^^v^  ^  y^-^^ 


=     -z 


ii  +  -21 


Z  9(^-1)  4^  -FM   , 
8y     X 


(14) 


If  ^I^V^  -Ic^^^^^ 


=   0  , 


(15) 


9c  -  e^'^  ' 


(16) 


3M 


9y     9c 


(17) 


where  M  =  uZ ,  M  =  vZ  and  H  =  0Z. 
X         y 

The  boundary  conditions  are 


C   =   0     at      c  =  0  , 


(18) 


and 


C   =   0     at      c  =  1  . 


(19) 
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For  this  study,  the  flow  will  be  independent  of  x  and  periodic 
in  y.   $   is  the  only  variable  allowed  to  vary  in  the  x- 
direction,  so  r —  =  0  for  all  other  variables. 

a  X 

In  order  to  obtain  an  expression  for  $,  equation  (16) 
may  be  integrated  with  respect  to  ^  giving 


$   =   /   -2.  e  z  dc  +  C  .  '  (20) 

0  ®0 


The  constant  of  integration,  C,  can  be  eliminated  by  taking 
the  vertical  average  of  equation  (20)  and  subtracting  from 
equation  (20)  which  gives  the  following 


$  =  -2_  J  e  zdc  -  /  e  zdc)  +  "^  .  (2i) 

^0   0  0 


In  the  above  equation,' ^represents  the  vertical  average  of 
a  given  variable. 

A  value  for  $  must  be  found  before  a  solution  for  $  can 
be  determined.   First,  take  the  vertical  average  of  (17)  and 
apply  the  boundary  conditions  (18)  and  (19)  which  gives 


I-  M^   =   -  Z  cl    =   0  .  (22) 


This  equation  shows  that  M  must  be  a  constant  in  order  to 
satisfy  the  boundary  conditions.   Next,  take  the  vertical 
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average  of  equation  (14)  and  differentiate  with  respect  to 
y  which  yields 


9   (-JT)  +   ^(^nP)       =   T-(-  z|?  +  ^(Ze)  (C-dH'  -  "f^J     .   (23) 


8t3y'  y'    ^^2'    y        8y^    8y    Bq'       '^^   '  9y      x 


The  first  term  in  equation  (23)  is  identically  zero  from  (22) 
and  the  resulting  equation  becomes 


h^^§     =  I?*-  Ij'V'  ^  e^'''"'^-"t-  ''"'  •   <'" 


This  equation  can  be  solved  diagnostically  for  $.   This  equa- 
tion must  be  applied  in  order  to  satisfy  the  boundary  condi- 
tions (18)  and  (19) . 

To  obtain  an  expression  for  -r— ,  which  is  a  constant, 

oX 

equation  (21)  is  differentiated  with  respect  to  x  which 
yields 


1^  =  ^^</o^^^-'/ll^^^'^f  • 


Since  -^  =  0  in  this  study, 

oX 


8_$   ^  ■d± 

ax    9x 


and  this  quantity  is  assumed  to  be  independent  of  y, 


(26) 
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Consider  the  domain  averages  of  equations  (13)  and  (14) 
which  are 


d  < — '        —  8'?     ' — ' 


vTl 


it<^'   =  -  ^  H  -  e^'^e)  (C-l)||  -  f  M^  ,  (28) 


8$ 
where  represents  the  y  average.   Initially,  :r—  is  computed 

geostrophically  as 


1^   =   ifT,  ,  .  (29) 

9x      7"   Y  t=0 


If  this  relationship  is  used  for  all  time  and  if  the  mass 

flux  in  the  x-direction  is  zero  initially,  then  M  will  be 

^  X 

zero  for  all  time.   The  first  two  terms  on  the  right  side 

of  equation  (28)  represent  the  mountain  drag.   The  sum  of 

these  two  terms  for  mid-latitude  flow  has  been  shown  to  be 

negative  in  the  mean.   This  will  lead  to  a  continual  loss  of 

y-momentum.   This  kind  of  loss  would  make  a  quasi-steady  final 

state  impossible.   An  artificial  source  of  momentum  is  required 

to  maintain  the  mean  meridional  flow.   This  artificial  source 

n 

of  momentum  is  provided  by  holding  tt —  constant  as  obtained 

oX 

in  equation  (29).   In  this  case,  if  "m   is  decreased  below 
the  initial  geostrophic  value,  M  will  increase  and  restore 
the  momentum  loss.   Thus,  in  the  long  term  mean,  the  mountain 
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terms  are  balanced  by  the  mean  Coriolis  force.   However,  for 
the  horizontal  scales  to  be  considered  in  this  study,  this 
momentum  source  may  not  be  needed  since  synoptic  scale  moun- 
tains cause  very  little  wave  drag. 

The  dependent  variables  are  to  be  divided  as  follows: 


u(x,y,C,t)   =   U  (x,y)  +  u(y,C/t)  , 


v(x,y,C/t)   =   V^(x,y)  +  v'(y,C/t)  , 

C(y/C/t)   =   C'(Y.C/t)  ,  (30) 

e(y,C/t)   =   e'(y,C/t)  , 

(|>(x,y,c.t)   =   $(x,y)  +  (j)'(y,C,t)  . 

The  basic  wind  deformation  field  is  chosen  to  be  periodic  in 
y  and  to  have  zero  vorticity  and  divergence.   Thus,  the 
deformation  field  is  given  by 

U^(x,y)   =   -D^  Lsinh(g)  sin  (^)  ,  (31) 


V,  (x,y)   =   -D,  Lcosh(^)  cos  (f )  ,  (32) 


where  D,  and  L  are  constants.   This  deformation  field  was 
d 

used  by  Stone  (1966)  in  his  study  of  quasi-geostrophic 
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frontogenesis.   Since  the  flow  is  assumed  to  be  independent 

of  X  in  this  study,  equations  (31)  and  (32)  are  evaluated  at 
X  =  0  resulting  in 

U^(0,y)   =   0  ,  (33) 


V^(0,y)   =-D^Lcos(^)  .  ^  (34) 


Substituting  the  expressions  in  equations  (30)  through  (34) 
into  equations  (13)  through  (15)  and  dropping  the  prime 
notation  yields 


8  { V  Z )       ^ — 

Z  1^  +  V^  z|^  +  u  —^—  +   l-(vZu)  +  |y(CuZ) 
dt     d   dy        dy      dy         9^ 


=   -  Z  1^  +  f  vZ  ,  (35) 

dX 


Z  1^  +  l-(vZv)  +  V  — ^  +  V.  z|^  +  It(^^2) 
9t    9y  dy       d   dy    dC 


=  -  z  1^  +  ^ez(c-i)  1^  -  fuZ  ,  (36) 

dy    Oq  dy 


^H^l^'-^e'  +!l<^9^>  =  °  '  <"> 


where  V^  =  -D ,  Lcos(^) .   This  last  set  of  equations  is  the 
d      d       L 

set  that  needs  to  be  solved  for  u,  v  and  6.   In  this  study, 
equation  (24)  does  not  have  to  be  solved  if  W'    is  restored 
to  the  initial  constant  value  at  every  timestep. 
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III.   FINITE  DIFFERENCE  FORMULATIONS 

Staggered  grids  are  used  in  both  the  horizontal  and 

vertical.   The  vertical  grid  is  divided  into  K-1  levels  of 

constant  A^  (Fig.  2).   Quantities  defined  for  the  layers  will 

be  denoted  by  integer  subscripts  and  the  ^  coordinate  value 

will  be  denoted  by  C^v*   The  quantities  9,  u,  v  and  cj)  are  all 

defined  on  the  integer  levels.   The  vertical  motion,  i, ,    on 

the  other  hand,  is  defined  on  the  top,  bottom  and  interface 

levels.   These  levels  will  be  denoted  with  half-integer 

subscripts  and  the  vertical  coordinate  will  be  identified  by 

^,  ,  .p.   The  horizontal  grid  is  divided  into  J-1  sections  of 

constant  resolution  (Fig.  3) .   Momentum  variables  are  defined 

on  the  integer  subscripted  points  denoted  by  Y . .   The  momentum 

variables  include  u,  v,  M  ,  M   and  V,.   The  mass  variables, 

X    y      d 

e,  (J),  Z  and  L, ,    are  defined  on  the  half-integer  subscripted 
points.   These  points  will  be  denoted  as  YI.  -,  ,2*   Values 
needed  off  their  respective  grid  points  will  be  identified 
by  a  caret  and  will  be  defined  as 


(Y.^,-Y.)S  .^T  /^  +  (Y.-Y.  -,)S.  T  /^ 
S    -     1+1   1   1+1/2      I 2^1 1-^/^  (38) 

n   ~  2 TYI      -YI       1  ' 

^  ^^^-^j+1/2  ^-^j-1/2^ 


for  quantities  defined  on  the  half-integer  points  and  as 
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^K1+1/ 


2        'C 


C\ 


K1 


»U,   V,   ^     <p 


«»K1  -V2 


^k  +  1 


k  +  1/2 


<l 


a    U,  V    6,    <P 


^k-1/2  "  C 


^3/2 

SI  /2 


a   U,  V,    ^,   <p 


Fig.  2.   The  vertical  grid,  showing  the  indexing 

convention  and  staggering  of  the  variables 
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S  -       '"i^3/2-^^.l/2'Si.l  ^  '"i^l/2-^Il-l/2'Si    ,3,^ 


for  quantities  at  the  integer  points.   Periodic  boundary 
conditions  are  used  in  this  study  for  all  quantities. 

Equations  (11),  (13)  through  (15)  and  (17)  are  solved  by 
using  finite  differences  in  y,  c,    and  t.   The  finite  differen- 
cing scheme  is  chosen  to  conserve  the  sum  of  kinetic  energy 
and  potential  energy. 

The  time  differencing  scheme  uses  a  leapfrog  scheme  with 
a  periodic  use  of  the  Euler  backward  scheme.   The  Euler 
backward  scheme  is  used  to  eliminate  the  computational  mode 
arising  with  the  leapfrog  scheme  [Haltiner  and  Williams, 
1980]  . 

Equations  (11),  (17)  and  (13)  through  (15)  may  be  written 
in  finite  difference  form  as 


MM-  MM 

^n+1  k"^  ^n  k  ^n+1  k"^^-;  k    ^i  k"^  ^i-l  k  ^i  k"^^n-l  k 

n          J  ,k         2  2                               2              2 

i     ^*  ^Ij+l/2  -  "j-1/2 


"i  k+l"^"-!  k      •  "-i  k'*"'\  k-1 

,  '^-^'j,k.l/2'^'   2  ^''^-  '^^'i,k-l/2'^'   2     ' 

'^k+1/2  ~  ^^k-1/2 


-^jH^^^jaV  ''°* 
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Z 


2j£  ^    '  2  '  '  2  '        ^  2  "  2 


J     9t  ^Vl/2-^Vl/2 


+   ^     j.k+1/2^  2  ^        ^^^j,k-l/2^  2  ^ 

^k+1/2    ~    ^k-1/2 


I     Vl/2A'^j-l/2.k        g    -  ^j+l/2-^j-l/2 


f-j,k^j,k'  ^^1) 


..  M^,  .^+3/2,k'^Qj+l/2,k.      M  Qj+l/2,k'^Qi-l/2,k. 


Z 


1+1/2  ,k         -j+l.k^ 2 ^        ^j.k 

j-M/2        3t  Vi""^: 


2 


^;7^                          .^i+l/2,k4-l'^^j+l/2,k,       •.                          ^j+l/2,k''^j+l/2,k-l 
_^    ^^^^j+l/2,k+l/2^  2  ^"^^^M+l/2.k-l/2^ 2 ^ 

^k+1/2  "   ^k-1/2 


=     0    ,  (42) 


^j+l/2,k+l         '^j  +  l/2,k      ^      _a.   7  ^j+l/2,k+l"^^j+l/2,k 

^k+l-^^k  "      ^0      ^^1/2  2 

(43) 

and 

My  -M 

j+l,k      yj,k       _    _  ^--^T  /-^    !'-■'  /-^    ^- 


/-j+l/2,k+l/2    ^j+l/2,k-l/2 
^■-     -    ^-  :  +  -L/^  '^k  +  1/2  ^k-1/2 


=    -Z 
j+1    -    "j  J'-/-  ^k  +  1/2  ''k-l/; 
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IV.   RESULTS 

All  numerical  results  to  be  shown  use  the  following 
numerical  values  for  the  constants : 

a   =   12  K 

B   =   375  m 

D   =   9  km 

D,   =   10~   sec"-^ 
a 

f   =   lO""^  sec"-"-  (45) 

2 
g   =   9.81  m/sec 

L   =   3  600  km 

e    =   300  K 


86 
s 

8z 


=   4  K/km 


The  following  finite-difference  increments  are  used. 

At   =   180  sec 

Ay   =   40  km  (46) 

Ac   =   1/27 

A.   MOUNTAIN  WITH  NO  DEFORMATION  FLOW 

The  initial  conditions  are  obtained  from  the  quasi- 
geostrophic  equations  for  steady  state  flow  over  topography 
[Walton,  1968;  DeBoer,  1970]  .   The  initial  potential 
temperature  field  is  given  by 
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e(y,z,o,  =  e^.B^^^iBMoU^^^g,  .2£_(J)  .  ,47, 

The  quantity  —r —  is  the  constant  initial  lapse  rate  and 

a  Z 

^  86^  1/2 


Substituting  equation  (4  7)  into  the  thermal  wind  equation 
yields 

^  =   1  -2-  i[B^  sinhla(z-D)]      y  _  2a      y  . 

8z      f  e^  L^^  3z    sinh(aD)    ^^"  ^L^    tt  ^^^  ^l^  ^  *       ^^^^ 

This  equation  can  be  integrated  with  respect  to  z  to  give 


/     r>\      T^  /  "^    s>  ^   cosh[a(z-D)]   .   -y, 
u  y,z,0)   =  B(^  -^— )     —  .'    )  sm  f 

■^  6^   dz        sinh(aD)        L 


-  ^  sin  (f)  (Z  -^)  .  (50) 

IT        L       z 


The  initial  v  field  can  be  shown  to  be 


/     A\      T7  ,   TTo  cosh[a(z-D)]     -y^  ,;--,. 

v(y,z,0)   =   V  +  aVB  — ^^^^j^^^^^y—  cos  (^)  ._       (51) 


The  initial  w  and  (p    are  obtained  from  equations  (4)  and  (3)  . 
It  is  necessary  to  rewrite  the  initial  conditions  in  the 
zeta-coordinate  system.   From  equation  (6),  equations  (47), 
(50),  and  (51)  can  be  rewritten  as 
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e„„,0,   =   e3.B!^iigIaZi^3,„,Z,  -2^3in(X,  ,  ,52) 


/     n^      i3^_g  ^^s,-'-/^  cosh[aZ(r-l)]   .   ,y, 
u(y,c,0)   =   ^^e^TT^       sinh(aD) ^^^  ^L^ 


-2^sin(^)(,-|)D  , 


(53) 


and 


v'y'C.O)   -   V  .  aVB  £2||ig2I^  COS  ,g,  .  (54) 

These  initial  conditions  are  shown  in  Figs.  4-7.   c  ^nd  (|)  are 
given  by  equations  (10)  and  (9). 

The  value  of  6   was  computed  using  the  formula 

^^S  D 

The  profile  of  the  mountain  used  in  this  study  is  given  by 


z    =   B  cos  (^)  .  (56) 

S  L 


and  is  shown  in  Fig.  8. 

The  initial  conditions  show  anticyclonic  vorticity  over 
the  ridge  and  cyclonic  vorticity  over  the  valley.  There  is 
divergence  on  the  upslope  and  convergence  on  the  downslope. 
These  perturbations  die  off  with  height  above  the  surface. 
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Figs.  9-12  show  the  results  at  t  =  72  hours.   Since  the 
original  fields  were  analytic,  only  a  slight  adjustment  to 
the  mountain  can  be  seen  in  the  v,  and  hence,  the  c  fields. 

B.   DEFORMATION  OVER  FLAT  TOPOGRAPHY 

For  this  case,  the  initial  potential  temperature  field 
contains  a  horizontal  variation.   The  initial  field  is  given 
by 


9  9 
e(y,z,0)   =  -y|-(z  -|)  -^cos(g)  .  (57) 


The  initial  u  field  is  again  derived  using  the  thermal  wind 
relationship  and  can  be  shown  to  be 


u(y,z,0)   =   -  ^  sin  (^)  (z  -|)  .  (58) 


The  initial  v  field  is  identically  zero  and  thus,  w  is  forced 
to  be  zero  by  equation  (3). 

It  is  again  necessary  to  write  equations  (48)  and  (49) 
in  zeta-coordinates  as 


9  S 
e(y,C.O)   =   -g^(c-l)D  -  ^cos  (^)  ,  (59) 


and 


u(y,C.O)   =   -  ^  sin  (g)  (c  -j)D  .  (60) 
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V  and  L,  are  zero  initially  and  (})  is  given  by  equation  (9). 
Figures  13  and  14  show  the  initial  conditions  and  Fig.  15 
shows  the  deformation  flow  given  by  equation  (34). 

The  front  is  developed  by  24  hours  as  shown  in  Figs.  16-19 
The  frontal  discontinuity  forms  along  the  point  where  v  =  0 
and  tilts  toward  the  colder  air.   Cyclonic  shear  can  be  seen 
along  the  frontal  boundary.   The  relative  vorticity  is  less 
than  zero  on  the  north  side  of  the  front  and  is  greater  than 
zero  to  the  south. 

By  48  hours.  Figs.  20-23  show  the  strengthening  of  the 
frontal  discontinuity.   The  tilt  remains  the  same  as  at  the 
24  hour  point.   On  the  other  side  of  the  domain,  the  tempera- 
ture gradient  is  weakening,  which  indicates  frontolysis 
occurring  in  that  region. 

C.   MOUNTAIN  WITH  NON-MOVING  DEFORMATION 

The  previous  two  cases  are  combined  in  order  to  investi- 
gate the  combined  effects  of  a  deformation  flow  and  a 
mountain.   The  initial  fields  are  given  by 


Q  /   ,-  A  \      A   _L  r»    s  sinh  [aZ  (^-1)  ]   .   ,y.   2a     ,y^    ,^-,, 
0(y,C,O)   =   9s  -^  ^  ^T    sinh(aD) ^^^  ^r^  "  T  ^°^  ^L^  '   ^^^^ 


r)  fi     1/2 

nfv  r  ni   -   Rr-2-  ^\  cosh[aZ(c-l)  ]  ^:,^,Y^ 

u(y,C,0)   -   B(Q^-y^)        sinh(aD) ^^^^L^ 


0 


-  %  sin  (f)  ,  (62) 

TT         Xj 
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v,y,c,0,      =      V  .   aVB   H|fi5t2|(piU  cos    ,J,      ,  ,63, 


and  ^  is  given  by  equation  (10) .   These  fields  can  be  seen 
in  Figs.  24,  25,  26  and  27.   The  deformation  is  the  same  as 
in  the  previous  case. 

As  in  the  previous  case,  a  tightening  of  the  temperature 
gradient  has  occurred  by  the  24  hour  point,  although  not  as 
strong  (Figs.  26-29).   The  front  has  moved  farther  to  the 
north  and  the  tilt  is  smaller  than  in  Case  B.   Throughout 
the  next  2  4  hours,  the  front  continues  to  move  to  the  north 
and  passes  over  the  mountain  ridge  (Figs.  30-33).   The 
gradient  weakens  as  the  front  passes  over  the  ridge  and  the 
tilt  is  negligible.   In  this  experiment,  the  front  weakens 
because  it  is  advected  away  from  the  frontogenetic  part  of 
the  deformation  field. 

D.   MOUNTAIN  WITH  MOVING  DEFORMATION 

The  next  four  cases  to  be  considered  start  with  the  same 
initial  conditions  as  in  Case  C  (Figs.  24,  25,  26  and  27). 
However,  in  these  cases,  the  deformation  moves  to  the  north 
at  a  speed  of  10  m/sec.   Each  case  examines  a  different 
starting  place  for  the  maximum  northward  deformation  flow. 
The  cases  have  the  maximum  flow  over  the  valley,  on  the 
upslope  side  of  the  mountain,  directly  over  the  ridge,  and 
on  the  downslope  side  of  the  mountain. 
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1.   Maximum  Deformation  over  the  Valley 

The  initial  deformation  is  shown  in  Fig.  15.  As  in 
the  previous  cases,  the  front  has  formed  by  24  hours  (Figs. 
34-37) .  There  is  a  small  tilt  toward  the  colder  air  in  the 
north  and  the  front  has  moved  farther  up  the  mountain  than 
in  Case  C.  The  upper  level  front  has  passed  over  the  moun- 
tain ridge. 

As  the  deformation  flow  continues  to  move  northward, 
so  does  the  front.   During  the  next  24  hours,  the  speed  of 
the  front  slows  down,  but  the  gradient  in  the  temperature 
field  tightens  (Figs.  38-41)  and  forms  a  discontinuity.   The 
cyclonic  shear  around  the  front  has  increased  and  there  is 
strong  vertical  motion  in  that  region. 

A  careful  comparison  of  the  potential  temperatures 
at  t  =  24  hours  between  Figs.  24  and  34  shows  very  little 
change  between  the  flat  experiment  and  the  mountain  experi- 
ment, except  for  the  translation  in  Fig.  34.   The  same  con- 
clusion holds  at  48  hours.   Bannon  (1983,1984)  has  shown  that 
mountain  effects  on  frontogenesis  will  be  important  when  the 
mountain  forced  convergence  is  of  the  same  order  as  the 
frontal  convergence.   Comparison  of  v  in  Figs.  6  and  18  shows 
that  the  mountain  forced  convergence  is  only  about  2  5  percent 
of  the  mountain  forced  convergence  at  24  hours.   This  ratio 
is  even  smaller  later  in  the  frontal  evolution.   With  the 
large  scale  mountain  used  in  this  and  later  experiments,  the 
effects  of  the  topography  are  unimportant. 
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2 .  Deformation  Shifted  One-Quarter  Wavelength 

In  this  experiment,  the  deformation  field  is  shifted 
one-quarter  of  a  wavelength  to  the  north  with  respect  to 
the  initial  fields  (Fig.  42) .   The  topography  is  also  unim- 
portant in  this  experiment.   At  t  =  24  hours,  there  is  only 
a  slight  indication  of  frontal  formation  in  the  temperature 
field,  but  a  much  stronger  suggestion  of  the  front  in  the  u 
field.   The  u  field,  however,  shows  little  tilt  (Figs.  43- 
46) .   Williams  (1968)  has  shown  that  this  phasing  of  the 
deformation  field  is  inefficient  for  frontogenesis. 

3.  Reversed  Deformation  Field 

As  is  shown  in  Fig.  47,  the  initial  deformation  flow 
is  in  opposite  phase  with  the  original  deformation  field. 
The  front  can  be  seen  at  24  hours  in  Figs.  48-51.   As  with 
previous  cases,  the  front  tilts  toward  the  cold  air,  but  in 
this  case,  the  tilt  is  toward  the  south.   This  is  a  result 
of  the  reversed  deformation  field.   The  fields  are  mirror 
images  of  the  fields  in  the  basic  experiment  except  for  the 
shift. 

4 .  Deformation  Reverse  of  Case  2 

The  last  case  to  be  considered  is  the  one  in  which 
the  maximum  deformation  flow  is  the  reverse  of  Case  2  (Fig. 
52) .   The  first  indication  of  frontal  activity  can  be  seen 
24  hours  later  (Figs.  53-56)  especially  in  the  u  field. 
Strong  cyclonic  activity  can  be  seen  here  with  no  evidence 
of  a  tilt.   The  temperature  field  shows  only  a  small  tightening 
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of  the  gradient  from  the  initial  conditions.   This  is  very 
similar  to  Case  2,  but  the  maximum  shear  is  at  the  lower 
boundary. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  study,  a  numerical  model  of  f rontogenesis  was 
constructed  which  included  surface  topography.   Periodic 
boundary  conditions  were  included  and  a  periodic  deformation 
field  was  introduced.   The  model  was  successfully  tested  with 
frontogenesis  with  no  topography  and  for  stratified  flow  over 
a  mountain  with  no  deformation  field.   Combined  experiments 
with  topography  and  a  moving  deformation  field  were  carried 
out.   The  experiments  showed  negligible  effects  on  fronto- 
genesis by  the  topography. 

The  next  step  in  this  research  is  to  introduce  smaller 
scale  topography  which  should  lead  to  a  noticeable  effect 
as  suggested  by  the  frontogenesis  theories  of  Bannon  (19  83, 
19  84)  .   Other  effects  including  vertical  wind  shear  could 
be  included  later. 
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